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ABSTRACT 

A recent theory is reviewed for the shot noise of coherent radiation propagating through 
a random medium. The Fano factor P j 1 (the ratio of the noise power and the mean 
transmitted current) is related to the scattering matrix of the medium. This is the optical 
analogue of Biittiker's formula for electronic shot noise. Scattering by itself has no effect 
on the Fano factor, which remains equal to 1 (as for a Poisson process). Absorption 
and amplification both increase the Fano factor above the Poisson value. For strong 
absorption P/I has the universal limit 1 + t|/, with / the Bose-Einstein function at 
the frequency of the incident radiation. This is the optical analogue of the one-third 
reduction factor of electronic shot noise in diffusive conductors. In the amplifying^ case 
the Fano factor diverges at the laser threshold, while the signal-to-noise ratio J 2 j P 
reaches a finite, universal limit. 

1. Introduction 

Analogies in the behavior of photons and electrons provide a continuing source of in- 
spiration in mesoscopic physics.u Two familiar examples are the analogies between 
weak localization of electrons and enhanced backscattering of light and between 
conductance fluctuations and optical speckle.cl The basis for these analogies is the 
similarity between the single-electron Schrodinger equation and the Helmholtz equa- 
tion. The Helmholtz equation is a classical wave equation, and indeed the study 
of mesoscopic phenomena for light has been limited mostly to classical optics. A 
common theme in these studies is the interplay of interference and multiple scatter- 
ing by disorder. The extension to quantum optics adds the interplay with vacuum 
fluctuations as a new ingredient. 

Recently a theoretical approach to the quantum optics of disordered media 
was proposed J3 that utilizes the methods of the random-matrix theory of quan- 
tum transport Jj'H The random matrix under consideration is the scattering matrix. 
The basic result of Ref. 3 is a relationship between the scattering matrix and the 
photocount distribution. It was applied there to the statistics of blackbody radia- 
tion and amplified spontaneoua^emission. This work was reviewed in Ref. 6. Here 
we review a later development^ the optical analogue of electronic shot noise. 



2 C. W. J. Beenakker and M. Patra 



Shot noise is the time-dependent fluctuation of the current I(t) = I + SI(t) 
(measured in units of particles/s) resulting from the discreteness of the particles. 
The noise power 



P = / dt6I(0)6I(t) (1) 

J — oo 

quantifies the size of the fluctuations. (The bar 7TT indicates an average over many 
measurements on the same system.) For independent particles the current fluctu- 
ations form a Poisson process, with power iioisson = I equal to the mean current. 
The ratio P/Pp isson (called the Fano factom) is a measure of the correlations be- 
tween the particles. 

For electrons, correlations resulting from the Pauli exclusion principle reduce P 
below Ppoisson- (See Ref. 9 for a review.) The ratio P/Pp issori-ig expressed in terms 
of traces of the transmission matrix t at the Fermi energy bytL3 

P/P -i (2 ) 

This formula holds at zero temperature (no thermal noise). In the absence of 
scattering all eigenvalues of the transmission-matrix product ffi are equal to unity, 
hence P = 0. This absence of shot noise is realized in a ballistic point contact JIIHI3 
At the other extreme, in a tunnel junction all transmission eigenvalues are <C 1, 
hence P = Pp i sson Jl3 A disordered metallic conductor is intermediate between these 



two extremes, having P = iPp i sson Jlil€ 

For the optical analogue we consider a monochromatic laser beam (frequency 
wo) incident in a single mode (labelled tuq) on a waveguide containing a disor- 
dered medium (at temperature T). The radiation from a laser is in a coherent 
state. The photostatistics of coherent radiation is that of a Poisson processJiS hence 
P = -Ppoisson for the incident beam. The question addressed in this work is: How 
does the ratio P/Pp isson change as the radiation propagates through the random 
medium? We saw that, for electrons, scattering increases this ratio. In contrast, in 
the optical analogue scattering by itself has no effect: P remains equal to Ppoisson 
if the incident beam is only partially transmitted — provided the scattering ma- 
trix remains unitary. A non-unitary scattering matrix, resulting from absorption 
or amplification of radiation by the medium, increases the ratio P/Pp isson- This 
excess noise can be understood as the beating of coherent radiation with vacuum 
fluctuations of the electromagnetic field.EZI 

Photon shot noise has been studied extensively in systemsj^here the scattering 
is one-dimensional (for example, randomly layered media) .I13't3 No formula of the 
generality of Eq. (||) was needed for those investigations. In order to go beyond the 
one-dimensional case, we have derived the optical analogue of Eq. (0). The result 



lSl 



P/Ppoisson = 1 + 2/(^0, r) [tt(1 it tti)t]m ° m \ (3) 

[t t\m m 

where f(ui, T) = [exp(hu> / kT) — I] -1 is the Bose- Einstein function. Eq. (S) contains 
both the transmission matrix t and the reflection matrix r (evaluated at frequency 



Photon Shot Noise 3 




Fig. 1. Coherent light (thick arrow) is incident on an absorbing or amplifying medium (shaded), 
embedded in a waveguide. The transmitted radiation is measured by a photodetector. 

ujq). For a unitary scattering matrix, rr^ + tft equals the unit matrix 1, hence 
the term proportional to / in Eq. (|J) vanishes and P — Pp \ sson . Absorption and 
amplification both lead to an enhancement of P above Ppoisson- For an absorbing 
system the matrix 1 — rr^ — W is positive definite and / > 0, so -P/Ppoisson > 1- 
In an amplifying system 1 — rr^ — tt' is negative definite but / is also negative 
(because T < in an amplifying system), so P/Pp isson is still > 1. 

We will review the derivation of the optical shot-noise formula (|^), and the 
application to absorbing and amplifying disordered waveguides. The ampli fy! tip 
case is of particular interest in view of the recent experiments on random lasers J 20 !' 2 H 
which are amplifying media in which the feedback required for a laser threshold is 
provided by scattering from disorder rather than by mirrors. 

2. Optical Shot-Noise Formula 

In this section we summarize the scattering formulation of the photodection problem,! 
and derive the formula (||) for the excess noise E We consider an absorbing or ampli- 
fying disordered medium embedded in a waveguide that supports N(u>) propagating 
modes at frequency cu (see Fig. [l]). The absorbing medium is in thermal equilib- 
rium at temperature T > 0. In the amplifying medium, the amplification could 
be due to stimulated emission by an inverted atomic population or to stimulated 
Raman scattering. til A negative temperature T < describes the degree of pop- 
ulation inversion in the first case or the density of the material excitation in the 
second case.ta A complete population inversion or vanishing density corresponds 
to the limit T — ► from below. The Bose-Einstein function f(u>,T) is > for 
T > and < —1 for T < 0. a The absorption or amplification rate l/r a = u\e"\ 
is obtained from the imaginary part e" of the (relative) dielectric constant (e" > 
for absorption, e" < for amplification). Disorder causes multiple scattering with 
rate l/r s and (transport) mean free path I — ct s (with c the velocity of light in 
the medium). The diffusion constant is D = ^cl. The absorption or amplification 
length is defined by £ a = \JDt & . 

The waveguide is illuminated from one end by monochromatic radiation (fre- 
quency to®, mean photocurrent Iq) in a coherent state. For simplicity, we assume 
that the illumination is in a single propagating mode (labelled mo). At the other 

a The quantity f(u>, T) is called the "population inversion factor" in the laser literature, because 
if u) is close to the laser frequency fi one can express / = (iVi owel ./JV U pp C r — I) 1 m terms of the 
ratio Wi ower /JVupper = exp(hQ/ kT) of the population of the lower and upper atomic levels, with 
/ = — 1 corresponding to a complete population inversion. 



4 C. W. J. Beenakker and M. Patra 



end of the waveguide, a photodetector detects the outcoming radiation. We assume, 
again for simplicity, that all N outgoing modes are detected with unit quantum ef- 
ficiency. We denote by p(n) the probability to count n photons within a time t. 
Its first two moments determine the mean photocurrent / and the noise power P, 
according to b 

I = -fi, P = lim - ( n 2 - n 2 ) . (4) 

The outgoing radiation in mode n is described by an annihilation operator 
a° ut (w), using the convention that modes 1,2, ... ,N are on the left-hand-side of 
the medium and modes N + 1, . . . , 2N are on the right-hand-side. The vector a out 
consists of the operators a° ut , a™*, . . . , a^w- Similarly, we define a vector a m for in- 
coming radiation. These two sets of operators each satisfy the bosonic commutation 
relations 

[a n (oj), at,(u/)] = S nm S(uj - u>'), [a„(w), a m (u>')} = 0, (5) 



and are related by the input-output relationst 

fl out = Sa in + Ub+VJ. (6) 

We have introduced the 2N x 2N scattering matrix S, the 2N x 2N matrices 
U, V, and the vectors b, c of 2N bosonic operators. The reflection and transmission 
matrices are N x TV submatrices of S, 

-(•;:')■ 

The operators b, c account for vacuum fluctuations. In order for these operators 
to satisfy the bosonic commutation relations (||) , it is necessary that 

Utf -VV ] = 1- SSl (8) 

In an absorbing medium c = and b has the expectation value 

(ftt (cj)b m (w')) = S nm 5(u - u')f(uo,T), T > 0. (9) 

Conversely, in an amplifying medium b = and c has the expectation value 

(c n (w)cj n (w , )) = -W(w-w')/(w,T), T<0. (10) 

The probability p(n) that n photons are counted in a time t is given byS 

p(n) = -(: W n e~ w :), (11) 



'This definition of P is equivalent to Eq. (|l|); In some papers the noise power is defined with an 
extra factor of 2. 
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where the colons denote normal ordering with respect to a out , and 

r t 2N 

W= dt> J2 a^tOCW. (12) 

J ° n=N+l 

a^\t) = {2-K)- 1 ' 2 due^a^H. (13) 
Jo 

Expectation values of a nor mal ly ordered expression are readily computed using the 
optical equivalence theorem.113 Application of this theorem to our problem consists 
in discretising the frequency in infinitesimally small steps of A (so that uo p = pA) 
and then replacing the annihilation operators a™(ui p ), b n (u> p ), c n (u> p ) by complex 
numbers a^p, b np , c np . The coherent state of the incident radiation corresponds to a 
non-fluctuating value of a™ p , such that |ajjp| 2 = S nmo S PPo 2ttIq / A (with uiq = poA). 
The thermal state of the vacuum fluctuations corresponds to uncorrelated Gaussian 
distributions of the real and imaginary parts of the numbers b np and c np , with zero 
mean and variance (\b np \ 2 ) = — (|c np | 2 ) = f(u> p , T), in accordance with Eqs. @ and 
©■ 

To evaluate the moments of the photocount distribution we need to perform 
Gaussian averages. The first two moments determine I and P. The results areSI 

p OO 7 

/ = Io[th] momo + — /(^^(l-rrt-ttt), (14) 



P = / + 27 /(wo ) 7 1 )[i t (l-rr t -tt^t] 



mama 



^/( W ,T) 2 Tr(l-rrt- tt t)2_ (15) 

The mean photocurrent is the sum of two terms, a term cx Io equal to the trans- 
mitted part of the incident current and a term cx / that represents the thermal 
emission of radiation. The noise power is the sum of three terms, the Poisson noise 
/ plus two sources of excess noise. The term cx f 2 is due to thermal emission while 
the term cx Iq/ is the excess noise due to beating of vacuum fluctuations with the 
incident radiation. For a unitary scattering matrix both terms vanish and P = I 
equals the Poisson value. 

The contributions from thermal emission to I and P can be eliminated by fil- 
tering the output through a narrow frequency window around ujq. Only the terms 
proportional to the incident current Iq remain, 

I = JbNmamo, (16) 

P = 7 + 2/ /(wo,T)[tt(l-rrt -tt^)t] momo . (17) 
This yields the optical shot noise formula (|^) discussed in the introduction. 

3. Absorbing Random Medium 

We consider an ensemble of absorbing disordered waveguides, with different realiza- 
tions of the disorder, and evaluate the ensemble averages of Eqs. ([l6]) and (|l^). For 
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Fig. 2. Excess noise power for an absorbing disordered waveguide, computed from Eqs. and 
Q). The ratio P/I tends to 1 + f / for L > £ a . 



a random medium the dependence on the index mo of the incident radiation is in- 
significant on average, so we may replace the average of a matrix element [• • -[mom,, 
by the average of the normalized trace N~ lr Tr. Moments of rr* and W in the 
presence of absorption have been computed by BrouwerEj using the methods of 
random-matrix theory, in the regime that both the length L of the waveguide and 
the absorption length £ a are much greater than the mean free path I but much less 
than the localization length Nl. This is the large- N regime N 3> L/l,£ a /l ^> 1. 
The ratio L/£ a = s is arbitrary. 
The result isS 



P 



Al s 
3L ° sinh s 



I 



21 J t 

M) IS 



2s + cotanh s s cotanh — 1 



sinh s 



sinh z 



sinh 



sinh 4 s 



(18) 
(19) 



The ratio P/Pp isson = P/I increases from 1 to 1 + |/ with increasing s, see 
Fig. ||. The limiting value P/Pp isson — > 1 + ^f(u>o,T) for L ^> ^ a depends on 
temperature and frequency through the Bose-Einstein function, but is independent 
of the scattering or absorption rates. This might be seen as the optical analogue 
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Fig. 3. Different dependence for electrons and photons of the shot-noise power P on the length 
L of the waveguide. The curve for photons is the same as in Fig. ^, the solid curve for electrons 
has been calculated in Ref. 25, the dashed curve is a qualitative interpolation. We have assumed 
a factor of 10 between the mean free path I, the absorption length £ a , and the localization length 
£l oc = Nl. (For electrons, the absorption length should be ignored.) The electronic P increases 
from to ^ of the Poisson value / when L becomes larger than I, and then increases further to 
full Poisson noise at £i oc = Nl. The photonic P has only a single transition, at £ a , from / to 
(1 + Nothing happens at L — I or L — £i oc to the shot noise of coherent radiation. 



of the_Lpiversal limiting value -P/Ppoisson — * g for L I of the electronic shot 
noise.00 

In the derivation of the limiting value of Pj I we have assumed that the length 
L of the waveguide remains small compared to the localization length £i oc = Nl. 
What happens in the localized regime L > £i oc ? Using the results from Ref. 24 we 
find that the ensemble averages (J) and (P) of current and noise in the localized 
regime are suppressed below the results Eqs. ( |T8| ) and ( |l9| ) in the diffusive regime, c 
but the ratio {P) / (I) remains equal to 1 + I/. This is a remarkable difference 
with the electronic analogue, where P becomes equal to the Poisson noise I in the 
localized regime. The difference between shot noise for electrons and photons is 
summarized in Fig. |^. 

c The precise result is {/) = (1 + |/) _1 (-P> = § (i/Ca)4 exp(-£/£ a - %L/NT) in the regime 
N, L/l > 5 a // > 1, for any value of L/Nl < 
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Fig. 4. Excess noise power for an amplifying disordered waveguide, computed from Eqs. ( pc| ) and 
). The ratio P/I diverges at the laser threshold L = 7r£ a . 



4. Amplifying Random Medium 

The results for an amplifying disordered waveguide in the large- N regime follow 



from Eqs. (|18[) and (19) for the absorbing case by the substitution r a 



-T a , or 



equivalenty s — > is. One finds 



I 

P 



4Z s 

or °~ 1 

AL sm s 



21 
ZL 



hfs 



2s — cotan s s cotan s — 1 



sm s 



■ 2 

sm s 



sin 3 s 



• 4 

sm s 



(20) 
(21) 



Recall that the Bose-Einstein function / < — 1 in an amplifying medium. As shown 
in Fig. H the ratio P/I increases without bound as the length L — > 7r£ a or, equiva- 
lently, the amplification rate l/r a — > ir 2 D/L 2 . This is the laser threshold. 

To understand better the behavior close to the laser threshold, we consider the 
scattering matrix S(u>) as a function of complex frequency uj. In the absence of 
amplification all poles (resonances) of S are in the lower half of the complex plane, 
as required by causality. Amplification shifts the poles upwards by an amount l/2r a . 
The laser threshold is reached when the first pole hits the real axis, say at resonance 
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frequency f2. For lo near f2 the scattering matrix has the generic form 

q cr n a m . . 

LO — \ l + — 1/^T a 

where <r n is the complex coupling constant of the resonance to the n-th mode in the 
waveguide and T is the decay rate. The laser threshold is at Tr a — 1. We will now 
show that, while P and I diverge at the laser threshold, Jjhe signal-to-noise ratio 
S = I 2 /P has a finite limit — independent of a n , T, or T a .Q 

We assume that the incident radiation has frequency u>q = f2. Substitution of 



Eq. fl22fl into Eqs. ( J20| ) and (21) gives the simple result 

2|/|E~ 



The total coupling constant E = Si + E r is the sum of the coupling constant 

£n=AH-l knP to * 

to 6 hence 



= En=i | cr T i 1 2 to the left end of the waveguide and the coupling constant E r = 
En=iv+i \ a n\ 2 to the right. The ensemble average (|ct„ 1o | 2 /S) is independent of 



since (Ej/S) = (£,/£) ^> (Ej/S) = 1/2. The signal-to-noise ratio of the incident 
coherent radiation (with noise power P = I ) is given by So = Iq/Pq = lo- The 
ratio S/Sq is the reciprocal of the noise figure of the amplifier. The signal-to-noise 
ratio of the transmitted radiation is maximal for complete population inversion, 
when |/| = 1 and (S) is smaller than So by a factor AN. This universal limit 

(S/So) — > 1/4JV does not require large AT, but holds for any N = 1, 2, It is the 

multi-mode gener a li s ation of a theorem for the minimal noise figure of a single-mode 
linear amplifier. I 



5. Outlook 

We conclude by mentioning some directions for future research. In the electronic 
case it is known that the result P/I = 1/3 for the Fajip factor of a diffusive conductor 
can be either computed from the scattering matrixL3 (using random-matrix theory) 
or from a kinetic equation known as the Boltzmann-Langevin equation.Ej Here we 
have shown using the former approach that the optical analogue is a Fano factor 
of 1 + |/ for a disordered waveguide longer than the absorption length. To obtain 
this result from a kinetic equation one needsa Boltzmann-Langevin equation for 
bosons. Work in this direction is in progress.tH 

The effect of localization on the Fano factor is strikingly different for electrons 
and photons. In the electronic case the average (P/I) goes to 1 in the localized 
regime, but we have found for the optical case that the ratio (P) / (I) of average noise 
and average current is unchanged as the length of the waveguide becomes longer 
than the localization length. We surmise that the same applies to the average (P/I) 
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of the ratio, but this remains to be verified. (In the diffusive regime the difference 
between the two averages can be neglected.) 

In the case of an amplifying disordered waveguide we have restricted ourselves 
to the linear regime below the laser threshold. Above threshold the fluctuations in 
the amplitude of the_£lectromagnetic field are strongly suppressed and only phase 
fluctuations remainES These determine the quantum-limited linewidth of the_radia- 
tion. A theory for this linewidth in a random medium is under development .E3 The 
application to a disordered waveguide would require a knowledge of the statistics 
of the poles of the scattering matrix in such a system, which is currently lacking. d 

The recent interest in the Hanbury-Brown and Twiss experiment for electrons in 
a disordered metaffl suggests a study of the optical case. The formalism presented 
here for auto-correlations of the photocurrent can be readily extended to cross- 
corrclationsJ23 but it has not yet been applied to a random medium. 

We do not know of any experiments on photon shot noise in a random medium, 
and hope that the theoretical predictions reviewed here will stimulate work in this 
direction. The universal limits of the Fano factor in the absorbing case and the 
signal-to-noise ratio in the amplifying case seem particularly promising for an ex- 
perimental study. 
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